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CHAPTER  7 

The  Relative  “Importance”  of 
Life-History  Stages  to  Population 
Growth:  Prospective  and 
Retrospective  Analyses 

Carol  Horvitz,  Douglas  W.  Schemske, 
and  Hal  Caswell 


Determining  the  importance  of  life-history  events  for  population 
growth  is  a  significant,  if  ill-defined,  goal  of  population-dynamics 
research.  Perturbation  analyses,  which  explore  the  effects  on  popu¬ 
lation  growth  of  changes  in  the  vital  rates,  provide  an  approach  to 
this  problem.  They  have  become  a  standard  part  of  demographic 
practice.  It  is  now  rare  to  find  a  published  report  of  population 
growth  rate  that  does  not  investigate  how  that  rate  changes  as 
the  vital  rates  are  perturbed,  either  actually  (comparing  differ¬ 
ent  treatments,  sites,  species,  etc.)  or  hypothetically  (exploring  the 
consequences  of  potential  management  strategies  or  of  evolutionary 
changes).  Applications  include  life-history  theory  (where  it  is  im¬ 
portant  to  know  how  the  different  vital  rates  influence  fitness;  see, 
e.g,,  Caswell  &:  Werner  1978;  Caswell  1985;  Calvo  &  Horvitz  1990; 
Kalisz  &:  McPeek  1992;  Calvo  1993),  conservation  biology  (where 
it  is  important  to  know  how  protecting  different  stages  in  the  life 
cycle  would  affect  population  growth;  see,  e.g.,  Crouse  et  al.  1987; 
Menges  1990;  Doak  et  al.  1994;  Heppell  et  al.  1994;  Schemske  et 
al.  1994),  ecotoxicology  (where  it  is  important  to  know  how  pollu¬ 
tants  affect  population  growth;  see,  e.g.,  Caswell  1996a;  Sibly  1996; 
Levin  et  al.,  in  press),  and  assessing  the  accuracy  of  estimates  of 
population  growth  rate  (Lande  1988). 
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Unfortunately,  as  the  use  of  demographic  perturbation  analyses 
has  spread,  so  have  some  common  misconceptions  about  them.  Our 
goal  in  this  chapter  is  to  clarify  some  of  the  applications  and  inter¬ 
pretations  of  perturbation  analyses.  Much  of  the  confusion  stems 
from  apparently  reasonable  but  ambiguous  questions,  perhaps  the 
most  ill  posed  of  which  is,  “which  of  the  stages,  or  which  of  the 
vital  rates,  is  most  important  to  population  growth?”  The  problem 
is  that  “importance”  has  several  different  meanings,  each  of  which 
leads  to  a  different  perspective  on  the  population. 

It  is  useful  to  distinguish  between  prospective  and  retrospective 
answers  to  this  question.  Prospective  analyses  address  the  effects 
of  potential  future  changes.  Of  all  the  changes  in  the  vital  rates, 
which  would  produce  the  biggest  effect  on  population  growth  rate? 
No  changes  in  any  vital  rates  need  have  occurred  to  ask  this  ques¬ 
tion.  It  is  even  possible,  and  sometimes  instructive,  to  ask  about 
changes  that  are  purely  hypothetical  (“if  pigs  had  wings  ...”  or  “if 
this  species  were  to  become  a  pereimial  instead  of  an  annual 
Sensitivity  and  elasticity  analyses  (see  Chapter  2  and  Caswell  1978; 
Caswell  et  al.  1984;  de  Kroon  et  al.  1986)  are  prospective  analyses. 

A  common  misconception  about  prospective  analyses  is  that  a 
high  sensitivity  or  elasticity  for  some  parameter  implies  that  an 
observed  change  in  population  growth  rate  was  due  to  that  param¬ 
eter.  This  is  a  mistake;  analyzing  changes  in  population  growth 
that  have  actually  occurred  is  a  retrospective,  not  a  prospective, 
problem.  Some  change  in  the  vital  rates  has  actually  occurred,  lead¬ 
ing  to  the  change  in  population  growth  rate.  We  want  to  know  how 
much  of  the  change  in  population  grovrth  rate  can  be  attributed 
to  the  changes  in  each  of  the  vital  rates.  The  vital  rate  that  con¬ 
tributes  the  most  to  the  variability  in  population  growth  is  not  nec¬ 
essarily  the  one  to  which  population  growth  rate  is  most  sensitive, 
nor  the  one  that  will  make  the  biggest  contribution  to  variability  in 
population  growth  rate  in  another  environment.  Some  vital  rates 
to  which  population  growth  is  very  sensitive  may  never  vary,  and 
thus  make  no  contribution  to  the  variability. 

The  appropriate  analytical  tool  for  this  question  is  the  “de¬ 
composition  analysis”  of  life-table-response  experiments  (Caswell 
1989a,  1996a).  If  population  growth  rate  is  measured  as  a  devia¬ 
tion  from  a  reference  value  (e.g.,  a  treatment  relative  to  a  control), 
then  the  treatment  effect  can  be  decomposed  into  contributions 
from  each  of  the  vital  rates.  If  variability  in  population  growth 
rate  is  expressed  as  a  variance  over  a  set  of  treatments,  this  vari¬ 
ance  can  be  decomposed  into  contributions  from  the  variances  and 
covariances  of  the  vital  rates. 
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Confusion  about  retrospective  and  prospective  perturbation  is 
analogous  to  some  common  misinterpretations  of  genetic-heritabil- 
ity  analysis  (Lewontin  1974).  Heritabiiity  decomposes  an  observed 
phenotypic  variance  into  genetic  and  environmental  contributions. 
It  is  a  retrospective  calculation.  High  heritabiiity  does  not  mean 
that  the  trait  is  insensitive  to  the  environment;  that  is  a  prospec- 
tive  conclusion.  The  trait  might  be,  for  example,  sensitive  to  tern- 
perature,  but  the  data  were  collected  in  a  constant-temperature 
chamber.  Similarly,  low  heritabiiity  does  not  mean  that  the  trait 
is  insensitive  to  changes  in  genetics.  Perhaps  the  trait  is  affected 
by  genotype,  but  the  population  studied  was  homozygous  at  the 
relevant  loci. 

We  describe  five  analyses,  two  prospective  and  three  retrospec¬ 
tive,  that  address  the  question  of  how  the  vital  rates  affect  popula¬ 
tion  growth  rate.  Then  we  apply  them  to  data  on  the  Neotropical 
understory  herb  Calathea  ovandensis.  Comparison  of  the  results 
sheds  light  on  various  aspects  of  the  “importance”  of  the  vital 
rates  in  this  particularly  well  studied  case. 


1  Demographic  Analyses 

We  use  the  linear  time-invariant  matrix  model 


n{t  +  1)  =  An{t) , 


(1) 


where  n{t)  is  a  vector  giving  the  abundances  of  the  stages  in  the 
population  at  time  t,  and  A  is  the  population-projection  matrix, 
whose  ijth  entry  a^j  gives  the  contribution  of  an  individual  in  stage 
j  to  stage  i  over  one  time  step.  See  Caswell  (Chapter  2,  or  1989a) 
for  details.  The  dominant  eigenvalue  A  of  the  matrix  gives  the 
population  growth  rate.  The  associated  right  and  left  eigenvectors 
w  and  V  give  the  stable  stage  distribution  and  the  stage-specific 
reproductive  values,  respectively. 

Sensitivity  and  elasticity  formulas  are  derived  in  Chapter  2.  The 
sensitivity  of  A  to  a  change  in  the  matrix  entry  aij  is  given  by 


_  9A  _  ViWj 

“  daij  ~  (w,v)  ’ 


(2) 


where  (w,  v)  is  the  scalar  product  of  w  and  v.  The  elasticity,  or 
proportional  sensitivity,  of  A  to  a  change  in  aij  is  given  by 


Ciij  ^A 

A  dcij  ' 


(3) 


It  is  sometimes  suggested  that  sensitivity  or  elasticity  analysis 
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is  better,  or  more  accurate,  or  less  biased,  or  more  biologically 
reasonable  than  the  other.  This  is  not  so.  The  two  give  different 
pictures  of  the  result  of  a  perturbation.  The  pictures  are  each  ac¬ 
curate,  unbiased,  and  biologically  meaningful.  Population  growth 
rate  (A)  is  a  function  of  the  entries  of  the  population-projection 
matrix.  Imagine  plotting  A  as  a  function  of  one  matrix,  say  a^-, 
while  holding  the  other  entries  constant.  Plot  aij  on  the  x-axis 
and  A  on  the  y-axis.  The  slope  of  the  resulting  curve  is  the  sen¬ 
sitivity  of  A  to  changes  in  aij.  If  you  were  to  repeat  this  opera¬ 
tion  for  all  the  matrix  entries,  the  result  would  be  a  set  of  curves 
with  different  slopes;  those  slopes  would  show  how  A  responds 
to  changes  in  all  the  matrix  entries.  Alternatively,  imagine  plot¬ 
ting  A  as  a  surface  in  a  multidimensional  space,  as  a  function 
of  all  the  aij.  The  sensitivities  Sij  give  the  gradient  of  this  sur¬ 
face,  showing  in  which  directions  it  slopes  steeply  and  in  which 
directions  it  is  relatively  fiat.  This  is  the  key  to  the  evolution¬ 
ary  applications  of  sensitivities,  which  can  be  interpreted  as  se¬ 
lection  gradients  (Lande  &  Arnold  1983;  Phillips  &  Arnold  1989). 

Now,  imagine  the  same  exercise,  but  plot  the  logarithm  of  A  as 
a  function  of  the  logarithm  of  aij.  The  slope  of  this  line  is  the 
elasticity  of  A  to  changes  in  aij.  The  multidimensional  analogue 
is  the  gradient  of  a  surface  that  shows  log  A  as  a  function  of  the 
logarithm  of  the  matrix  entries.  The  elasticity  Cij  =  0  if  aij  =  0. 
That  is,  if  the  wingspan  of  a  pig  increases  by,  say,  10  percent, 
there  is  no  change  in  its  flying  ability.  Transitions  directly  from 
newborn  to  large  reproductive  individuals  (entries  in  the  lower  left 
corner  of  the  matrix)  usually  have  high  sensitivities,  even  though 
such  transitions  may  not  occur.  The  sensitivity  indicates  correctly 
what  would  happen  if  that  rate  were  increased  from  zero  to  a  small 
positive  number.  The  elasticity,  however,  is  zero  and  can  reveal 
nothing  about  such  changes.  This  makes  elasticity  inappropriate 
for  evolutionary  questions  such  as  asking  how  fitness  would  change 
if  an  annual  were  to  become  perennial,  or  seed  dormancy  were 
to  be  introduced  where  it  does  not  exist,  etc.  The  elasticity,  in 
contrast,  indicates  correctly  that  increasing  the  matrix  entry  by 
some  proportion  has  no  effect  on  A. 

The  elasticities  sum  to  1  across  the  whole  matrix  (Caswell  1986; 
de  Kroon  et  al.  1986;  Mesterton-Gibbons  1993)  and  can  be  in¬ 
terpreted  as  proportional  contributions  of  the  corresponding  vital 
rates  to  A  (see  van  Groenendael  et  al.  1994).  There  is  some  confu- 
^sion  between  the  correct  statement  “stage  x  contributes  y  percent 
to  population  growth  rate”  and  the  incorrect  statement  “stage  x 
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explains  y  percent  of  the  variation  in  population  growth  rate,” 
which  cannot  be  determined  from  an  elasticity  analysis. 

Is  elasticity  more  meaningful  than  sensitivity?  The  choice  of 
arithmetic  or  logarithmic  axes  is  a  matter  not  of  right  or  wrong  but 
of  which  graph  reveals  the  aspects  of  the  curve  that  are  of  interest. 
Logarithmic  axes  are  useful  because  equal  intervals  correspond  to 
equal  proportions.  If  you  are  interested  in  patterns  involving  pro¬ 
portional  changes  in  vital  rates,  use  the  slopes  of  the  logarithmic 
plots  (i.e.,  the  elasticities).  Equal  intervals  on  arithmetic  axes  im¬ 
ply  equal  changes.  If  you  are  interested  in  rates  of  change  on  a 
linear  scale  (as  evolutionary  questions  are,  for  example),  then  use 
sensitivities.  Better  yet,  make  both  calculations  and  compare  the 
results- 

Remember  that  both  sensitivity  and  elasticity  are  derivatives; 
they  give  the  local  slope  of  A  (or  log  A)  as  a  function  of  aij  (or 
log Oij).  They  may  not  accurately  predict  the  result  of  large  per¬ 
turbations  (although,  in  practice,  they  do  a  surprisingly  good  job). 
Using  the  second  derivatives  (Caswell  1996&)  to  take  into  account 
the  curvature  of  A  might  improve  the  prediction.  A  more  straight¬ 
forward  approach,  however,  is  to  vary  parameters  and  calculate  A 
numerically  (see,  e.g.,  Caswell  &  Werner  1978;  Bierzychudek  1982; 
Horvitz  &  Schemske  1986;  Martinez  &  Alvarez-Bullya  1986;  Calvo 
1993).  For  example,  an  effective  but  rare  pollinator  may  increase 
plant  reproductive  success  several-fold.  A  local-perturbation  anal¬ 
ysis  near  the  usual  low  level  of  fruit  production  may  not  accurately 
predict  the  effect  of  the  pollinator  on  A,  while  a  simulation  that 
actually  varies  fruit  set  and  calculates  A  will  do  so  (Horvitz  & 
Schemske,  unpubl.  data).  Numerical  perturbation  analysis  is  also 
useful  for  investigating  specific,  biologically  interpretable  changes 
in  several  of  the  vital  rates  simultaneously.  It  should,  however,  be 
supplemented  with  analytical  results. 


2  Retrospective  Analysis 

A  retrospective  analysis  begins  with  data  on  the  vital  rates  and 
on  A,  under  tw’o  or  more  sets  of  environmental  conditions.  The 
goal  of  the  analysis  is  to  quantify  the  contribution  of  each  of  the 
vital  rates  to  the  variability  in  A.  A  formal  method  for  quantifying 
data  of  this  kind  has  been  called  Life-Table-Response  Experiments 
(LTRE’s)  (Caswell  1989a,  6,  1996a,  in  press;  see  also  Levin  et  al. 
1987,  in  press;  Walls  et  al.  1991;  Brault  k  Caswell  1993).  This 
usage  defines  “experiment”  in  a  sense  wide  enough  to  include  not 
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only  designed  manipulations  but  also  comparative  observations. 

Intuitively,  variation  in  one  matrix  entry,  say  a^j,  makes  a  small 
contribution  to  the  variation  in  A  if  aij  doesn’t  change  much,  or 
if  A  is  not  very  sensitive  to  a^j,  or  if  both  are  true.  Thus,  the 
contribution  of  aij  to  the  variation  in  A  involves  the  product  of  Sij 
and  the  observed  variation  in  aij. 

As  in  the  familiar  analysis  of  variance,  LTRE  analysis  depends 
on  whether  the  “treatments”  are  considered  a  random  sample  from 
some  universe  of  possible  conditions  or  a  fixed  set  of  treatments 
that  are  of  interest  in  themselves  and  could  be  repeated  if  nec¬ 
essary.  As  in  statistics,  the  distinction  between  random  and  fixed 
treatments  is  subtle,  and  one  data  set  may  well  be  interpreted  dif¬ 
ferently  by  different  people.  In  a  random  design,  the  results  are 
characterized  by  the  variance  of  A,  and  the  goal  of  the  analysis 
is  to  decompose  this  variance  into  contributions  from  the  vari¬ 
ances  of  (and  covariances  among)  the  matrix  entries  (Brault  & 
Caswell  1993;  Caswell  &  P.  Dixon,  unpubl.).  In  a  fixed  design,  the 
results  are  described  in  terms  of  the  effect  of  each  treatment  on 
A,  measured  relative  to  some  baseline.  The  goal  of  the  analysis 
is  to  decompose  the  treatment  effects  into  contributions  from  the 
treatment  effects  on  each  of  the  vital  rates.  (There  are  also  LTRE 
methods  for  quantitative  treatments,  analogous  to  regression  mod¬ 
els;  see  Caswell  1996a;  Caswell  &  L.V.  Martin,  unpubl.) 


Random  Treatments:  Variance  Decomposition 

Let  V  (A)  denote  the  variance  of  A  among  treatments.  Recall  the 
formula  for  the  variance  of  a  linear  combination  of  two  random 
variables  x  and  y: 

V{ax  -\-by)=  a^V{x)  +  b'^V{y)  +  2ahC{x,  y),  (4) 

where  C{x,y)  is  the  covariance  of  x  and  y.  Approximating  A  as 
a  linear  function  of  the  aij,  and  using  the  sensitivities  Sij  as  the 
slopes,  yields 

VW  > ki)sijske,  (5) 

ij  ki 

where  C(ij,  ki)  is  the  covariance  of  aij  and  and  the  sensitivities 
are  calculated  at  the  mean  matrix.  The  covariances  are  calculated 
directly  from  the  data  on  A  for  each  environment.  Each  term  in 
this  summation  is  the  contribution  of  one  pair  of  vital  rates  to 
y(A).  This  calculation  was  first  used  to  examine  the  variation  in  A 
among  pods  of  killer  whales  (Brault  &  Caswell  1993). 
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Because  y(A)  depends  on  the  covariances  among  pairs  of  vital 
rates,  it  is  difficult  to  say  which  rate,  as  opposed  to  which  pair 
of  rates,  makes  the  biggest  contribution.  One  way  to  define  the 
contribution  of  a  single  rate  is  to  sum  the  contributions  for  all  the 
covariances  involving  that  rate: 

Xij  =  ^  C{ij,  ktjSijSu-  (6) 

H 

This  sum  includes  the  contribution  to  1^(A)  from  the  variance  of 
Oij  plus  half  the  contributions  from  the  covariances  of  with 
the  other  rates.  The  other  half  of  the  covariance  contributions  is 
distributed  among  the  contributions  of  the  other  matrix  entries. 
We  call  this  calculation  of  the  contribution  Xij  tfis  “covariance 
method." 

A  special  case:  independent  variation.  If  it  seems  safe  to  assume 
that  the  vital  rates  vary  independently,  then  the  approximation  for 
the  variance  reduces  to 

V{X)^J2Via,j)sl,  (7) 

where  V{aij)  is  the  observed  variance  in  the  ijth  element  of  the 
matrix  A.  Each  term  in  this  summation  represents  a  contribution 
of  one  vital  rate  to  y(A);  by  definition, 

^  =  (8) 

As  with  the  covariances  in  the  preceding  case,  the  variances  here 
are  calculated  from  the  observed  set  of  matrices  and  the  sensitivi¬ 
ties  are  calculated  at  the  mean  matrix.  We  call  this  the  “variance 
method”  of  calculating  the  contribution  Xij  to  V{X). 

Choice  of  a  reference  matrix.  The  variance  V (A)  is  the  average 
of  the  squared  deviations  from  the  mean  of  the  growth  rates,  and 
thus  the  sensitivities  in  the  approximations  (6)  and  (8)  are  calcu¬ 
lated  from  the  mean  matrix.  However,  other  matrices  might  serve 
as  more-reasonable  reference  matrices.  In  the  data  set  we  analyze 
next,  we  use  a  summary  matrix  computed  from  the  pooled  data 
from  all  sites  and  years.  This  is  equivalent  to  using  a  weighted  mean 
with  weights  proportional  to  the  sample  sizes.  In  this  case,  V(A) 
is  no  longer  strictly  speaking  a  variance:  it  is  the  mean-squared 
deviation  from  the  reference  matrix,  and  is  no  longer  a 

covariance;  it  is  the  mean  of  the  cross-product  of  the  deviations  of 
Oij  and  Oki  from  their  values  in  the  reference  matrix. 
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Fixed  Treatments:  Decomposing  Treatment  Effects 

A  fixed-effect  analysis  treats  the  matrices  as  representative  of  par¬ 
ticular  conditions,  either  experimental  or  natural  (high  vs.  low  nu¬ 
trients  in  a  one-way  model,  for  example,  or  year  and  spatial  loca¬ 
tion  in  a  two-way  model).  The  goal  is  to  determine  how  much  of 
the  main  effect  of  each  treatment  level  on  A  is  contributed  by  each 
of  the  vital  rates.  In  a  factorial  design,  this  can  be  extended  to 
include  the  interaction  effects  of  combinations  of  factors  as  well  as 
the  main  effects. 

The  analysis  uses  a  linear  approximation  in  which  the  sensitiv¬ 
ities  appear  as  slopes.  The  effect  of  a  treatment  on  A  depends  on 
its  effect  on  each  matrix  entry  and  on  the  sensitivity  of  A  to  that 
entry. 

In  order  to  analyze  a  data  set  in  which  plots  and  years  appear 
as  fixed  effects,  we  write  a  two-way  model  (as  in  an  analysis  of 
variance)  for  population  growth  in  plot  m  and  year  n: 

^(mn)  _  ^(..)  _|_  ^  (0'/3)rnn  ,  (9) 

where  A^"^  is  the  growth  rate  calculated  from  ,  the  grand  mean 

of  all  the  matrices,  and  /?n,  and  {a(3)mn  are  the  plot,  year,  and 
interaction  effects.  These  effects  are  estimated  by 

-  A(-) ,  (10) 

/?„  =  A(")-AW,  (11) 

(aPUr.  =  (12) 

Here,  the  growth  rates  and  A^'^^  are  the  growth  rates  calcu¬ 
lated  from  the  mean  over  years  of  the  matrices  for  plot  m  and  from 
the  mean  over  plots  of  the  matrices  for  year  n,  respectively. 

These  treatment  effects  on  A  can  be  decomposed  into  contribu¬ 
tions  from  the  effects  on  each  matrix  element,  as  follows: 

Oi„.  «  ,  (13) 

ij 

,  (14) 

ij 

(aP)mn  «  -  a|v^)sij  '  (15) 

ij 

(Caswell  1989a, 6).  Because  the  sensitivity  structure  may  change 
from  one  treatment  to  another,  the  sensitivity  matrix  is  evaluated 
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at  a  matrix  that  is  intermediate  between  the  treatment  being  con¬ 
sidered  and  its  reference  matrix.  For  ami  this  intermediate  matrix 
is  +  A^‘))/2;  for  it  is  (A(-^)  +  A^-))/2;  for  (a/?)^n,  it 

is  (A(^’")+A("))/2. 

As  in  the  random-effects  case,  it  is  possible  to  use  matrices  other 
than  the  mean  as  a  reference  matrix.  In  our  example,  instead  of  the 
overall  mean,  we  use  the  summary  matrix  from  the  pooled  data  for 
A^"\  To  be  consistent,  we  should  also  have  used  pooled  data  for  the 
plot  means  A^^'^  and  the  year  means  A^”^^  We  did  not  have  the 
data  accessible  in  this  format  at  the  time  of  these  analyses,  so  we 
did  not  do  this.  Thus,  we  present  the  analysis  with  the  cautionary 
note  that  our  reference  matrices  axe  slightly  inconsistent,  but  we 
are  reasonably  confident  that  this  does  not  have  a  large  effect  on 
the  qualitative  results. 

The  decomposition  analysis  (eqs.  13-15)  gives  the  contribution 
of  each  matrix  entry  to  each  treatment  effect  (i.e.,  to  each  plot 
effect,  to  each  year  effect,  and  to  each  plot-by-year  interaction  ef¬ 
fect).  It  would  be  nice  to  have  some  way  to  compute  a  single  num¬ 
ber  that  measures  the  overall  contribution  of  each  matrix  entry 
to  each  treatment  factor.  Some  matrix  entries  are  unaffected  by 
the  treatment  and  make  a  zero  contribution  to  each  level  of  the 
treatment.  Others  make  small  positive  contributions  at  some  lev¬ 
els  and  small  negative  contributions  at  other  levels.  The  important 
matrix  entries  are  those  with  large  positive  contributions  at  some 
treatment  levels  and  large  negative  contributions  at  others.  Tak¬ 
ing  the  means  of  these  contributions  over  the  different  levels  of  the 
factor  does  not  work,  because  it  is  always  approximately  zero,  just 
as  in  an  analysis  of  variance  the  mean  of  the  treatment  effects  is 
(exactly)  zero. 

This  suggests  using  the  variance,  or  perhaps  the  mean  of  the 
absolute  values,  of  the  contributions  as  a  summary  measure.  In 
this  chapter  we  use  the  means  of  the  absolute  values.  This  cal¬ 
culation  is  useful  in  simplifying  the  results  of  complex  data  sets 
with  many  levels  of  cross-classified  factors.  In  our  system,  with 
4  plots,  4  years,  and  16  plot-year  combinations  (see  below),  tak¬ 
ing  these  means  yields  a  graph  with  a  single  line  for  the  plot  ef¬ 
fects,  a  single  line  for  the  year  effects,  and  a  single  line  for  the 
interaction  effects.  However,  this  simplification  comes  at  a  cost: 
presenting  the  results  in  this  form  obscures  the  contributions  to 
the  treatment  and  interaction  effects,  which  are  the  original  goal 
of  the  design  (for  complete  factorial  analyses,  see  Caswell  1989a; 
Walls  et  al.  1991). 


256 


Horvitz,  Schemske  &  Caswell 


Table  1.  Variance  ofX  among  16  Matrices  for  Calathea 
ovandensis 


Approximation 

Equation 

Value 

Variance  from  mean 

standard 

0.0182 

Variance  from  pooled 

16 

0.0187 

Covariance  approximation 

5 

0.0317 

Variance  approximation 

7 

0.0196 

3  An  Example:  Calathea  ovandensis 

We  now  apply  the  analyses  described  above  to  a  set  of  16  matri¬ 
ces  (four  plots  over  four  years)  generated  from  a  five-year  study 
of  a  Neotropical  xmderstory  herb,  Calathea  ovandensis.  The  life 
cycle  was  divided  into  eight  stages:  seeds,  three  vegetative  stages 
(seedling,  juvenile,  and  pre-reproductive)  and  small,  medium,  large, 
and  extra-large  reproductive  plants.  C.  ovandensis  has  long-term 
seed  dormancy.  Seedlings  have  the  highest  mortality  of  all  stages, 
and  reproductives  have  low  mortality.  Plants  that  do  not  die  some¬ 
times  shrink,  sometimes  remain  the  same  size  (stasis),  and  some¬ 
times  increase  in  size  between  seasons.  Fertility  is  positively  corre¬ 
lated  with  size  (Horvitz  &  Schemske  1995). 

Significant  spatiotemporal  variability  is  found  in  the  vital  rates 
and  in  the  sensitivity  structure  (Horvitz  &  Schemske  1995).  The 
population  growth  rates  of  the  16  plot-year  matrices  varied  from 
0.7356  to  1.2477,  with  an  arithmetic  mean  of  0.9695  and  a  variance 
y(A)  =  0.0182  (Table  1).  Nevertheless,  the  overall  dynamics  are 
well  represented  by  a  summary’  matrix  that  pooled  all  obser¬ 
vations  over  the  five-year  period  (Fig.  1).  The  dominant  eigenvalue 
of  this  matrix  is  =  0,9923  (Horvitz  &  Schemske  1995). 

We  characterize  the  variability  in  the  population  growth  rate  by 
the  mean-squared  deviation  around  the  growth  rate  of  the  sum¬ 
mary  matrix,  rather  than  as  the  variance  around  the  mean.  Let 
A(Tnn)  dgjjote  the  growth  rate  in  plot  m  in  year  n.  Then,  our  mea¬ 
sure  of  variability  is 

14(A)  =  ^(A^’^^)  -  A^-))Vl6 . 


(16) 
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Figure  1.  Entries  aij  of  the  pooled  summary  matrix  (Horvitz  & 
Schemske  1995)  for  Calathea  ovandensis.  The  64  entries  of  the  S  x  S 
matrix  are  arranged  in  column  order:  fates  of  seeds  fan, . . .  ,a8i)j  fates 
of  seedlings  fdi2, . . .  ^as2),  and  so  on.  Subscripts  for  the  entries  of  the 
first  row  appear  on  the  abscissa;  subscripts  for  the  diagonal  elements 
appear  directly  above  the  matrix  entry. 

The  variance  from  the  summary  matrix  calculated  in  this  way  is 
y,(A)  =  0,0187  (Table  1). 


Prospective  Analyses 

The  sensitivities  and  elasticities  of  the  summary  matrix  are  shown 
in  Figure  2.  The  largest  sensitivities  are  sgi,  571,  sei,  and  S51,  in 
that  order.  Additive  perturbations  thus  exert  the  biggest  effects  if 
they  occur  in  the  vital  rates  that  represent  extremely  rapid  growth 
of  newly  germinated  seedlings.  Sensitivities  to  other  seed  transi¬ 
tions  (S41  and  S31)  and  to  the  rapid  growth  of  established  seedlings 
(^82,  572,  562,  552  and  542)  are  also  high.  Within  each  stage,  stasis 
and  growth  have  higher  sensitivities  than  fecundity. 

The  largest  elasticity  is  655  (stasis  of  small  reproductives),  fol¬ 
lowed  by  654  (growth  of  pre-reproductives),  644  (stasis  of  pre-repro- 
ductives),  and  en  (seed  dormancy).  Thus,  proportional  perturba¬ 
tions  have  the  biggest  impact  if  they  affect  the  stasis  of,  and  growth 
to,  small  reproductives. 

There  is  little  correlation  between  the  sensitivity  and  the  elas¬ 
ticity  patterns  (Fig.  2).  This  is  not  unusual;  the  effects  of  additive 
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Matrix  parameter  subscript,  ij 


Figure  2.  (a)  Sensitivities  (sij)  ofX  to  the  corresponding  entries  of  the 
pooled  summary  matrix,  (b)  Elasticities  (dj )  of  \  to  the  corresponding 
entries  of  the  pooled  summary  matrix.  Axis  labels  as  in  Figure  1. 

and  proportional  perturbations  are  likely  to  differ. 

Retrospective  Analyses:  Random  Effects 

The  covariances  of  the  matrix  entries  are  shown  in  Figure  3a.  The 
large  positive  values  dominating  the  figure  are  covariances  among 
the  fecundities;  a  year  or  plot  that  is  good  for  the  reproduction 
of  one  size  class  tends  to  be  good  for  the  reproduction  of  all.  The 
entries  on  the  diagonal  of  this  surface  are  the  variances;  these  are 
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Figure  3.  (a)  Covariances  of  the  matrix  elements  aij  and  aki  over  16 
plot-year  combinations.  Variances  appear  on  the  diagonal,  (b)  The  prod¬ 
ucts  Si j Ski  of  pairs  of  sensitivities^  calculated  from  the  pooled  summary 
matrix,  (c)  Contributions  of  each  pair  of  matrix  entries  to  Vs  (A),  from 
the  product  of  the  corresponding  elements  in  (a)  and  (b) .  The  order  of 
the  matrix  elements  in  all  three  surfaces  is  the  same  as  the  abscissa  in 
Figure  1. 
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Figure  4.  The  variances  of  the  matrix  elements  Oij  over  16  plot-year 
combinations.  Axis  labels  as  in  Figure  1. 

plotted  in  Figure  4.  The  variances  of  fecundities  dwarf  all  the  other 
variances. 

The  products  SijSu  of  pairs  of  sensitivities  are  highest  for  the 
germination  and  rapid  grovi:h  of  seeds  (Fig.  36).  Thus,  in  this  data 
set  there  is  a  negative  correlation  between  the  sensitivity  of  A  to 
a  matrix  entry  and  the  variability  in  that  matrix  entry.  The  con¬ 
tributions  to  T4(A)  from  each  of  the  covariances  is  given  by  the 
element-wise  product  of  Figures  3a,  6;  the  result  is  shown  in  Figure 
3c.  Several  large  contributions  j&rom  variance  terms  appear  on  the 
diagonal,  and  there  are  many  small  but  not  negligible  contribu¬ 
tions,  both  positive  and  negative,  from  the  off-diagonal  covariance 
terms.  The  sum  of  these  values,  as  in  equation  (5),  is  0.0317. 

Taking  sums  across  the  rows  (or  down  the  columns)  of  the  surface 
in  Figure  3c  gives  the  contributions  Xij  by  the  covariance  method 
(eq.  6).  The  results  are  shown  in  Figure  5.  This  shows  how  (A) 
is  increased  by  positive  covariances  (including  variances)  of  matrix 
entries  and  reduced  by  negative  covariances  among  matrix  entries. 
The  largest  of  these  contributions  are  from  (growth  of  small  re- 
productives),  as4  (growth  of  pre-reproductives),  031  (rapid  gro'wth 
of  newly  germinated  seeds),  and  a45  (a  negative  contribution  from 
covariances  involving  the  shrinking  of  small  reproductives).  Nega¬ 
tive  contributions  are  also  made  by  covariances  involving  a44  and 
a55  (stasis  of  pre-reproductives  and  small  reproductives). 

Figure  3  shows  that  covariances  among  the  matrix  entries  can¬ 
not  be  ignored  in  this  data  set.  Nevertheless,  Figure  56  shows  the 
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(a) 


(b) 


Figure  5.  (a)  The  total  contribution  (xij)  of  matrix  entry  aij  to  the 
variance  14(A),  calculated  by  summing  over  the  covariance  contributions 
in  Figure  3  according  to  equation  (6).  (b)  The  contributions  (Xij^) 
the  variance  of  matrix  entries  to  14(A)  assuming  that  the  matrix  entries 
vary  independently j  as  in  equation  (8).  Axis  labels  as  in  Figure  i. 


contributions  calculated  using  the  variance  method  (eq.  8); 
these  values  are  just  the  diagonal  elements  of  Figure  3c.  They  sum 
to  0.0196  (Table  1). 

The  largest  contributions  to  14  (^)  analysis  are  from  the 

variances  of  055  (growth  of  small  reproduct ives),  asi  (growth  of 
newly  germinated  seeds),  a55  (stasis  of  small  reproductives),  and 
^154  (growth  of  small  reproductives).  Note  that  ignoring  covariances 
leads  to  the  conclusion  that  the  variance  of  ass  makes  a  large  pos- 
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Year 


Figure  6.  Population  growth  A  for  each  plot  in  each  year,  expressed  as 
a  deviation  from  the  growth  rate  of  the  pooled  summary  matrix. 

itive  contribution  to  14(A);  in  fact,  the  covariances  involving 
make  its  net  contribution  to  14(A)  negative  (Fig.  56). 


Retrospective  Analyses:  Fixed  Treatments 

Using  the  fixed-effect  factorial  design,  we  can  break  down  contribu¬ 
tions  to  variability  in  A  further,  into  independent  effects  of  plots, 
years,  and  their  interaction.  Population  growth  rate  is  shown  in 
Figure  6  as  a  function  of  plot  and  year.  The  plot  and  year  effects 
am  and  are  shown  in  Figures  7a, 6.  Plot  2  had  the  largest  posi¬ 
tive  effect  on  A,  and  plot  4  the  largest  negative  effect.  Year  1982-83 
had  the  largest  positive  effect,  and  1984-85  the  largest  negative  ef¬ 
fect.  The  plot  and  year  effects  are  of  similar  magnitude,  although 
years  are  slightly  more  influential  than  plots.  The  temporal  pattern 
differed  among  plots  (Fig.  6),  suggesting  a  plot-by-year  interaction, 
which  appears  in  the  interaction  effects  {al3)mn  in  Figure  7c. 

We  do  not  show  the  decomposition  of  each  of  the  plot,  year, 
and  interaction  effects  into  contributions  from  each  of  the  matrix 
entries.  Instead,  we  show  the  means  of  the  absolute  values  of  the 
contributions  of  each  matrix  entry  to  the  plot  effects,  to  the  year 
effects,  and  to  the  interaction  effects  (Fig.  8). 

The  matrix  entries  most  influential  in  determining  plot  effects 
are  aes  (growth  of  small  reproducti ves) ,  044  (stasis  of  pre-reproduc- 
tives),  031  (germination  and  growth  of  seeds),  and  045  (shrinkage 
of  small  reproductives)  (Fig.  8a). 

The  entries  most  influential  in  determining  year  effects  axe  055 


IMPORTANCE  OF  LIFE-HISTORY  STAGES 


263 


Figure  7.  (a)  The  main  effects  am,  m ’=■  oj  plots,  calculated 

from  equation  (13).  (b)  The  main  effects  (Sn,  n  =  of  years, 

calculated  from  equation  (U)-  (c)  The  interaction  effects  {al3)mn,  cal¬ 
culated  from  equation  (15).  Years,  from  1982  to  1985,  are  shown  within 
each  plot. 
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(a) 


Figure  8.  (a)  The  mean,  over  4  plots,  of  the  absolute  value  of  the  con¬ 
tributions  of  each  matrix  element  to  the  plot  effect  (b)  The  mean,  over 
4  years,  of  the  absolute  value  of  the  contributions  of  each  matrix  ele¬ 
ment  to  the  year  effect  (c)  The  mean,  over  16  plot-year  combinations, 
of  the  absolute  values  of  the  contributions  of  each  matrix  element  to  the 
plot-by-year  interaction  effect  Axis  labels  as  in  Figure  1. 


IMPORTANCE  OF  LIFE-HISTORY  STAGES 


265 


(stasis  of  small  reproductives),  031  (germination  and  growth  of 
seeds),  045  (regression  of  small  reproductives),  and  a2i  (seed  ger¬ 
mination)  (Fig.  86). 

The  matrix  entries  with  the  largest  influence  on  interaction  ef¬ 
fects  are  045  (growth  of  pre-reproductives),  aes  (growth  of  small 
reproductives),  031  (germination  and  growth  of  seeds),  and  ass 
(stasis  of  small  reproductives)  (Fig.  8c). 


4  Discussion 

We  show  computations  of  seven  different  measures  of  the  “im¬ 
portance”  of  individual  matrix  entries  to  population  growth:  two 
prospective  and  five  retrospective.  Here  we  discuss  the  matrix  ele¬ 
ments  that  rank  among  the  top  four  in  each  type  of  analysis  (Table 
2).  We  also  examine  the  correlations  among  the  seven  measures  of 
importance  (Table  3). 

The  sensitivity  and  elasticity  analyses  pick  out  completely  differ¬ 
ent  sets  of  matrix  entries  as  the  most  important.  Additive  pertur¬ 
bations  of  the  matrix  elements  in  the  last  four  locations  of  the  first 
column,  representing  extremely  rapid  growth  of  newly  germinated 
seedlings,  have  the  biggest  effect  on  A.  Proportional  perturbations, 
in  contrast,  have  the  biggest  impact  when  they  involve  the  sta¬ 
sis  and  growth  of  pre-reproductive  and  small  reproductive  plants, 
or  seed  dormancy.  Sensitivity  is  weakly,  but  negatively,  correlated 
with  indices  generated  by  the  other  approaches  (Table  3),  in  spite  of 
the  fact  that  sensitivity  appears  in  the  calculation  of  all  the  other 
indices.  This  correlation  would  change  if  it  was  calculated  using 
only  the  sensitivities  of  the  nonzero  matrix  entries,  that  is,  those 
that  are  free  to  vary.  The  largest  such  sensitivity  is  531,  and  we 
note  that  asi  also  appears  among  the  four  most  important  entries 
in  all  five  retrospective  analyses. 

Of  the  64  matrix  entries,  only  eight  (ass,  ^^54,  ciu,  0315 

a4s,  and  a2i)  appear  among  the  top  four  in  importance  according 
to  elasticity  or  any  of  the  retrospective  analyses.  None  is  important 
by  all  criteria.  The  growth  of  newly  germinated  seedlings  (asi)  is 
important  by  all  five  retrospective  indices;  it  contributes  signifi¬ 
cantly  to  variance  (according  to  both  the  covariance  and  variance 
methods),  and  its  contribution  to  plot,  year,  and  plot-by-year  in¬ 
teraction  effects  is  also  important. 

Only  two  matrix  entries  (ass  and  aes,  the  stasis  and  growth  of 
small  reproductives)  are  important  according  to  four  indices.  Entry 
ass  has  high  elasticity  and  makes  a  large  contribution  to  T4(A)  by 
the  variance  method.  It  has  an  important  contribution  to  the  year 
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effect  and  to  the  plot-by-year  interaction  effect  (Table  2).  Entry 
a65  makes  a  large  contribution  to  V"5(A)  and  to  the  plot  and  plot- 
by-year  interaction  effects. 

Shrinkage  of  small  reproductives  (045)  and  the  growth  of  pre- 
reproductives  (a54)  are  important  according  to  three  of  the  indices. 
Shrinkage  of  small  reproductives  has  high  elasticity  and  makes  im¬ 
portant  contributions  to  Vs  (A)  by  the  variance  method  and  to  the 
plot-by-year  interaction.  The  growth  of  pre-reproductives  makes  an 
important  contribution  to  Vs{X)  and  to  the  plot  and  year  effects. 

Stasis  of  pre-reproductives  (044)  has  high  elasticity  and  makes 
an  important  contribution  to  the  plot  effect.  Seed  dormancy  and 
seed  germination  (an  and  a2i)  have  high  elasticity  and  make  an 
important  contribution  to  the  year  effect,  respectively. 

Table  3  shows  the  correlations  between  the  various  measures 
of  “importance,”  underscoring  the  different  insights  produced  by 
each  analysis.  Sensitivity  and  elasticity  are  nearly  independent  of 
contributions  to  V5(A)  calculated  using  equations  (5)  and  (6),  the 
appropriate  variance  decomposition  for  this  data  set  (correlation 
coefficients  —0.07  and  0.01,  respectively).  Thus,  treating  the  16 
matrices  as  a  random  sample  characterized  by  the  variance  of  A, 
the  prospective  and  retrospective  analyses  clearly  address  different 
questions.  Elasticity,  however,  is  positively  correlated  with  contri¬ 
butions  to  Vs{X)  by  the  variance  method. 

Sensitivity  is  negatively  correlated  with  the  results  of  all  the 
other  analyses.  Elasticity  has  positive  correlations  with  most  other 
analyses,  and  large  ones  with  the  plot  (0.64),  year  (0.77),  and  plot- 
by-year  interaction  (0.69)  effects. 

The  random-  and  fixed-effects  analyses  are  positively  correlated, 
but  not  always  strongly.  Contributions  to  Vs{X)  from  the  covariance 
approach  has  correlations  of  0.41,  0.14,  and  0.47  with  contributions 
to  the  plot,  year,  and  plot-by-year  interaction  effects.  Within  the 
fixed-effect  design,  there  are  strong  correlations  between  contribu¬ 
tions  to  plot  and  year  effects  (0.72),  plot  and  interaction  effects 
(0.89),  and  year  and  interaction  effects  (0.74).  This  suggests  that 
the  plot  and  year  effects  tended  to  be  mediated  by  variation  in  the 
same  set  of  vital  rates  (Table  2). 

Finally,  we  note  that  is  highly  correlated  with  elasticity  and 
with  contributions  to  all  three  fixed  effects.  Thus,  a  random-effects, 
retrospective  analysis  that  is  inappropriate  for  this  data  set  does  a 
good  job  of  predicting  the  results  of  both  a  fixed-effect,  factorial, 
analysis  and  an  elasticity  analysis. 

There  is  no  single  best  method  for  determining  the  “importance” 


It 


IMPORTANCE  OF  LIFE-HISTORY  STAGES 


269 


of  the  vital  rates.  It  is  important  to  define  more  precisely  what  we 
mean  by  importance,  and  to  address  specific,  well-posed  questions, 
such  as  the  following. 

What  are  the  relative  sensitivities  of  A  to  incremental  changes 
in  the  vital  rates? 

WRat  are  the  relative  proportional  sensitivities  of  A  to  propor¬ 
tional  changes  in  the  vital  rates? 

Over  an  observed  sample  of  environmental  conditions,  how  do 
the  variances  of,  and  covariances  among,  the  vital  rates  con¬ 
tribute  to  the  variance  of  A? 

How  much  does  each  of  the  vital  rates  contribute  to  the  observed 
effect  on  A  of  each  of  a  set  of  fixed  treatments? 

In  the  case  analyzed  here,  the  prospective  elasticity  analysis 
identifies  one  vital  rate,  seed  dormancy,  that  is  not  important  ac¬ 
cording  to  the  retrospective  methods.  Similarly,  the  factorial  anal¬ 
ysis  identifies  seed  germination  as  making  an  important  contribu¬ 
tion  to  year  effects;  this  vital  rate  is  not  identified  as  important  by 
any  other  method.  Despite  these  differences,  both  the  transition 
to  and  the  fates  of  small  reproductives  are  identified  as  key  events 
for  Calathea  ovandensis,  both  in  predicting  the  outcome  of  hypo¬ 
thetical  perturbations  and  in  explaining  the  observed  variation  in 
population  growth  rate. 
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